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EQUATION  FOR  USE  IN  INTERNAL  BALLISTIC  CALCULATIONS 


Prepared  by: 
A*  E.  Selgel 


ABSTRACT:  The  solution  of  the  internal  ballistics  problem  requires  a 

knowledge  of  the  equation  of  state  of  the  propellant  gases.  Usually 
tho  Noble -Abel  equation  6f  state  is  used  for  this  purpose;  However, 
this  equation  does  not  account  for  the  mtersolecular  forces. 

A seed-empirical  equation  of  state  has  been  developed  which 
permits  the  solution  of  the  internal  ballistics  problem  and  repre- 
sents accurately  the  effects  of  the  intersolecular  forces  over  a wide 
range  of  pressures  and  densities-.  Since  the  semi  -empirical  equation 
is  similar  in  fora  to  the  perfect  gas  equation,  it  is  as  convenient 
in  application  to  interior  ballistic  calculations. 
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This  report  presents  some  of  the  results  of  a theoretical  and  experi- 
mental study  made  at  the  University  of  Amsterdam  and  the  Nrval  Ordnance 
Laboratory  with  an  experimental  gun,  the  "Expansion  Rate  Measuring 
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A CONVENIENT  AND  ACCURATE  SEMI-EMPIRICAL  ENTROPIC 
EQUATION  FOR  USE  IN  INTERNAL  BALLISTIC  CALCULATIONS 


I , INTRODUCTION 

In  order  to  calculate  the  mu2zle  velocity  of  a gun,  the  iaentropic 
equation  of  state  of  the  powder  gases  must  be  known.  For  this  purpose 
the  ballisticians  universally  employ  the  Nobie-Abel  iaentropic 
relation , 


p(v-b)  = constant, 

where  b and  8 are  assumed  constant  during  the  expansion  and  are 
the  co-volume  and  mean  ratio  of  specific  heats  of  the  propellant, 
respectively.  This  equation  is  based  upon  ths  Abel  thermal  equation 
of  state,  p(v-b)  = RT,  and  the  assumption  that  the  specific  heat  is 
constant.  Although  it  is  an  approximate  representation  of  the 
behavior  of  propellant  gases,  its  use  is  justified  qy  the  fact  that 
the  approximation  is  aa  good  as  others  introduced  in  the  theory  of 
internal  ballistics  and  that  it  is  mathematically  extremely  tractable. 
Nevertheless,  it  would  be  desirable  to  have  an  isentropic  equation 
which  more  accurately  describes  ths  behavior  of  the  propellant  gases 
and  which  is  convenient  for  ballistic  calculations. 

Such  an  en tropic  relation  for  a real  gas  (one  in  whloh  the  inter- 
molecular  forces  are  important,  in  contrast  to  a perfect  or  Abel  gas 
in  which  ths  in termol ocular  forces  are  unimportant)  was  found  by 
the  present  author  in  his  study  of  the  rapid  expansion  of  compressed 
gases  behind  s piston.  (This  study  was  made  at  the  van  der  Uaals 
Laboratories,  Amsterdam,  and  is  described  fully  in  reference  (a ).) 

Ths  entropic  relation,  which  is  semi-empirical,  and  the  manner  in 
which  it  was  obtained  are  discussed  in  this  paper. 


II.  THE  EXPERIMENTAL  GUN  APPARATUS 

Ths  experimental  gun  apparatus  employed  in  this  work  was  first 
used  by  Michels,  Slawsky,  and  Jacobs  (reference  (b));  it  Is  desig- 
nated aa  ERMA  (from  the  initial  letter?  of  the  descriptive  title 
"Expansion  Rate  Measuring  Apparatus").  The  instrument  is  essentially 
a hollow  cylinder,  closed  at  one  end  and  open  at  the  other  end. 

Inside  the  constant  cross-sectional  area  cylinder  is  a movable  piston 
or  projectile  of  known  macs  which  is  initially  restrained  from 
movement.  Between  the  piston  snti  tho  closed  end  is  contained  a 
highly  compressed  gas  whose  temperature  and  pressure  are  known. 

During  the  operation  of  the  experiment  the  compressed  gas  propels 
the  piston  from  rest  out  of  the  op* n snd  of  the  cylinder.  From  the 
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experiment  the  travel  and  time  measurements  of  the  piston  could  be 
determined.  These  crave! -time  data  were  used  to  study  the  dynamic 
expansion  process  in  the  experimental  gun. 

ill.  THE  BALLISTIC  THFORY 

The  basic  assumptions  used  to  describe  the  motion  in  ERMA  of  the 
gas  and  the  projectile  are  the  following: 

(1)  The  gas  motion  is  one- dimensional  in  space. 

(2)  Each  part  of  the  gas  expands  isantropically. 

(3)  The  projectile  is  propelled  by  the  gas  pressure  and  opposed 
by  a frictional  resistance  between  the  projectile  and  the  cylinder 
walls  and  by  the  air  pressure  in  front  of  the  projectile. 

Proa  these  assumptions,  hydrodynamic  and  thermodynamic  lavs  lead 
to  the  "characteristic  equations".  It  can  be  then  shown  that  two 
types  of  expansions  are  possible.  One  type  is  the  "non-reflection* 
case  lit  whioh  the  behavior  of  the  piston  is  the  same  as  would  occur 
if  the  back  wall  of  the  chamber  were  infinitely  distant  from  the 
piston;  the  second  type  ia  the  "reflection"  case  in  which  the  behavior 
of  the  piston  is  influenced  by  the  presence  of  the  back  end  of  the 
chamber.  For  both  cases  the  entire  behavior  of  the  expending  gas  and 
moving  piston  can  be  determined  from  the  characteristic  equations, 
the  real  gas  equation  of  state  (in  the  form  of  an  entropic  pressure- 
density  relationship),  and  from  Newton's  foroe  lav  applied  to  the 
projectile. 


IV.  OBTAINING  THE  SEMI-EMPIRICAL  ENTROPIC  EQUATION 

As  a result  of  tests  performed  vitn  the  experimental  gun  on 
nitrogen  gas,  it  was  previously  observed  that  the  use  <">f  "iaplifisd 
equations  of  state,  such  as  the  perfect  gas  or  ..oci*s  equating, 

yielded  unsatisfactory  agreement  between  the  oellist'  ineorv  and  cue 
experiment;  satisfactory  agreement  was  only  obtained  »>•  isen- 

tropic  nitrogen  aata  «a  determined  fres  experiment-'  ; ;w-T  measurements 
was  used  (see  references  (a)  and  (b)). 

An  attempt  was  made  to  obtain  an  analytic  expression  relating  the 
pressure  behind  the  projectile  with  the  projectile  velocity  for  the 
uon-refleotion  (i.e.,  infinite  chamber)  case.  Therefore , from  the 
p- data  of  nitrogen,  Ison  tropic  curves  of  piston  velocity,  u. , 
as  a function  of  pressure,  p , were  obtained  for  the  non-reflection 
case.  Proa  these  curve  a derived  only  from  nitrogen  thermodynamic 
date,  the  position- time  relation  of  the  projectile  could  be  obtained 
in  the  iaentropic  case  of  unopposed  projectile  from  the  expressions, 
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Various  relations  were  attempted  to  approximate  the  p-u-  curves 
(power  series,  rational  fractions,  etc.)  without  success.  Here  the 
difficulty  lay  in  obtaining  an  expression  which  fitted  the  curves  and 
yet  could  be  conveniently  integrated  to  yield,  for  example,  a relation 
between  projectile  velocity  and  barrel  length. 

It  was  recalled  that  the  relation  between  p and  UL  for  a perfect 
gas  in  the  non-reflection  case  is 


3t«  A * - 0 

X- 1 _a  \ 

2 <W 


where  p0  is  the  initial  pressure  of  the  expansion,  ouQ  is  the  initial 
sound  velocity,  and  V is  the  ratio  of  specific  heats.  In  view  or 
this  a similar  expression  for  the  real  gas  was  attempted. 


where  ft  is  assumed  a function  of  entropy  alone  and  c<  0 is  an  initial 
parameter  of  the  expansion.  Then  from  the  condition  that  equation  (2) 
must  hold  with  a constant  f>  along  an  isentropic  p - a curve  for  all 
initial  preasures  p0,  follows  the  fact  that  is  equal  to  <J>  pj* 

Efrttt  $ W&rgET-gigBS* 

This  semi-empirical  equation  satisfies  the  requirements  of  a 
desired  analytic  expression  very  well.  With  the  proper  values  of  £ 
and  4>  it  fits  the  isentropic  p-LL curves  within  the  accuracy  with  which 
these  curves  were  obtained  (i  part  in  50  to  1 part  in  500),  and  it  is 
extras^  ly  convenient,  sines  it  readily  integrated  as  necessary. 


/ros  equations  (1)  and  (2)  ttz  expression  for  projectile  travel  as 
a function  of  time  for  an  unopposed  projectile  propelled  by  an  isen- 
troplcaily  expanding  gas  in  the  non-reflection  case  car.  be  obtained: 

f P -2 

x - *.t  - -rrrfnfe —P.t.  i-  il'  — l>o) 
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are  X = C 
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Equation  (2)  c-*n  h*  written  an  a relation  between  the  Rieoann 
function  cr  (defined  below)  and  the  gaa  pressure  along  an  isentropc. 


P = P„(l  + ) 


* 


(2a) 


where  the  Riemann  function  O" 
equation  (2a)  ia  independent  of 
relationship  between  gas  properties. 


It  is  to  be  noted  that 


perimsntal  gun;  it  expresses  a 


By  use  of  the  definition  of  C the  above  expression  can  be  trans- 
formed to  obtain  a relation  between  pressure  and  density  for  i sen tropes. 
The  following  result  is  obtained. 


g-‘Vr  i , 
P If 


P(P~2) 


where  f is  a function  of  entropy  alone  (as  is  K).  It  is  seen  that 
f is  the  volume  occupied  by  a gas  if  compressed  i sen tropically  to  an 
infinite  pressure.  Equation  (4)  is  the  ssmi-eapirical  antropic  equation. 


V.  METHOD  OF  OBTAINING  THE  CONSTANTS  OF  THE  EQUATION 

Each  entropy  has  associated  with  it  different  valucc  of  p , ^ . and  f ; 
these  are  dependent  only  on  the  entropy  and  the  gas.  To  obtain  these 
constants  for  a given  entropy,  equation  (4)  was  used  with  p-  9 isen- 
tropic  data.  At  three  points  on  the  lsentropic  p - 9 curve  corresponding 
values  of  p and  Q were  obtained.  Then  by  equation  (4) 


and 
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which  was  solved  for  p ; $ was  then  obtained  from  equation 
from  equation  (4).  A plot  ia  given  in  figure  1 of  p , 4>  , 
as  a function  of  entropy  for  nitrogen  gas. 


(5),  and  f 
and  f 
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VI.  ACCURACY  AND  RANGE  OF  THE  EQUATION 

The  range  In  which  the  experimental  gun  expansion  occurs  extends 
free;  states  at  about  50°C  and  high  pressures  (400  atmospheres  to  300C 
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■ii  to  states  at  low  tear  •-stares,  -?5WC  to  -12  '°C,  ana 
lev  pressures  (of  about  100  atmospheres;  ':i..  isentropes  S = -12 
through  S = -17  eai/aole°C  for  nitrogen  cover  this  range  ( S is  chosen 
to  be  equal  to  zero  at  1 atmosphere  and  0°C).  Therefore,  the  constants 
for  the  seod-empirical  equation  were  calculated  for  these  entropy 
values. 


The  average  difference  between  densities  calculated  from  the  semi- 
empirical  equation  (4)  and  those  read  fro*  a plot  of  nitrogen  isentropes* 
is  0.25£  down  to  200  atmospheres.  Between  200  atmospheres  and  100 
atmospheres  the  deviation  is  greater,  becoming  1%  in  some  cases  at  100 
atmospheres.  It  is  to  be  noted  that  in  this  low  pressure  region  the 
accuracy  of  the  data  is  also  less  (about  l/2£) , and  the  co-existence 
region  is  beginning  to  influence  the  isentropes. 

Although  the  constants  of  the  semi -empirical  equation  for  nitrogen 
were  chosen  to  apply  to  regions  below  50°C,  the  equation  has  been 
employed  successfully  (1%  deviation  in  density)  at  125°0  and  600C 
atmospheres  with  these  same  constants.  Thus,  it  is  probable  that  the 
equation  can  be  applied  to  higher  pressure  and  temperature  regions. 


For  the  only  i sen trope  examined  which  extended  into  a perfect 
g*s  region,  the  smai-empirical  equation  seemed  to  apply.  It  fits 
the  nitrogen  isentrope,  5 = -4  cal/mole°C,  in  the  range  of  -125°C 
and  1 atmosphere  to  +150°C  and  34  atmospheres.  The  equation  has  been 
applied  with  success  to  argon  in  the  region  of  0°C  to  150°C  and  pres- 
sures up  to  2000  atmospheres. 


It  seems  likely  that  the  semi -empirical  equation  could  be  applied 
to  powder  gaa  in  ballistic  problems  with  success. 


VII.  THERMODYNAMIC  RESULTS  OF  THE  EQUATION 


From  the  en tropic  equation  of  state,  equation  (4) , should  follow 
all  the  thermodynamic  quantities.  Some  of  these  are  listed  below: 


P^(v  - 
U=  q 

H = 9 


x X 

' " P(jl-E) 

+ «Z/2P 


X. 


+ fp  + 


2 

c* 


2(p-2) 


en tropic  equation  of  state 


irternal  energy  equation 


enthalpy  equation 


*Theae  isentropes  were  taken  from  reference  (c^,  where  the  thermo- 
dynamic properties  of  nitrogen  for  pressures  up  to  6000  atmospheres 
between  -125°C  and  +150°C  are  given. 
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£ P 

o<p 


velocity  of  sound 


V 4 

. d 3 ^ df  temperature  equation 

dS  + Pd5 

g , f , 8 , 4>  t 11,10  are  functions  only  of  entropy,  and 
jeciflc  volume  = 1/f  , CX  = 4>p*'^  • 


VIII.  THE  USE  OF  THE  SEMI -EMPIRICAL  EQUATION  TO  OBTAIN 
THE  DESCRIPTION  OF  THE  EXPERIMENTAL  GUN  PERFORMANCE 

If  the  semi-empirical  equation  is  used,  a relation  between  pro- 
jectile relocity  and  travel  can  be  obtained  for  an  isentropic  expansion 
when  there  are  no  reflections  which  reach  the  projectile. 


P„A  v _ I 

j[p-S)(p-0 

Thus,  the  diitensionless  projectile  relocity  U/o<o  can  be  plotted 
against  the  dimensionless  projectile  position  p0Ay./M  o'.*  to  yield 


fp-e)  ~ fe -»)(<-£ 

(•  - *.r 


curves  for  each  entropy  of  performance  for  the  isentropic  non- 
reflection unopposed  projectile  case.  These  curves  are  shown  as  solid 
lines  in  figure  2 for  nitrogen  gas. 

If  reflections  reach  the  projectile,  these  curves  of  course  are 
altered.  By  using  the  semi-empirical  equation  and  the  numerical  char- 
acteristics method  the  dimonslonless  velocity  can  be  obtained  as  a 
function  of  dimensionless  projectile  position  for  a given  G/M 
(where  G>  Is  the  mass  of  gas  behind  the  piston,  M is  the  projectile 
mass)  and  entropy.  In  figure  2 these  G/M  lines  are  plotted 
dashed,  to  indicate  the  reflected  cases. 

It  was  found  from  tests  made  with  the  experimental  gun  that  the 
forces  opposing  the  projectile  could  be  conveniently  and  accurately 
accounted  for  by  the  use  in  the  calculations  of  an  "effective*  pro- 
jectile mass,  M/0  —T^) , in  place  of  the  actual  projectile  mass,  M. 

(see  references  (a)  end  (d)h  Thus,  the  curves  of  figure  2 can  be 
applied  to  the  opposed  projectile  case  (if  everywhere  the  projectile 
case  M is  replaced  by  M/(l-*f\,)  ) and  are  a complete  description  of 
the  experimental  gun  performance. 
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IX.  APPLICATION  OF  THE  Sa<I-B*PIAICAL.  EQUATION  TO 
ACTUAL  BALLISTIC  CALCULATIONS 

The  seal -empiri cal  equation  can  be  as  conveni  sr.tly  applied  to 
ballistic  calculations  as  the  Noble-Abel  or  perfect  gas  eauations. 

The  necessary  gas  constants  for  such  calculation,  6 and  $ , should 
be  obtained  by  fitting  the  seal -empirical  isentropic  equation  to  the 
powder  gas  p~9  iaentropic  data.  If  this  data  is  unavailable,  p 
and  4>  can  be  approximately  found  from  the  data  of  tvo  gun  firings 
made  at  different  values  of  dimensionless  length  or  velocity  by  the 
use  of  the  non-reflection  isentropic  theory,  i.e.,  from  equations  (3) 
and  (6).  (The  application  of  the  non-reflection  isentropic  theory  to 
gun  firings  baa  been  found  to  be  a viseful  approximation,  especially 
to  high  velocity  gun  firings • See  reference  (e ) . ) 
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LIST  OF  SYMBOLS 

GL  - Velocity  of  sound 

A - Cross-sectional  area  of  projectile 

brtf  - Co-volume  in  Abel  thermal  equation  of  state 

r 

f - Entropy  constant  in  semi-empirical  en tropic  equation 

'0  - Entropy  constant  in  internal  energy  and  enthalpy  equations 

G - Mass  of  gas  behind  the  projectile 

H - Enthalpy 

M - Mass  of  projectile 

p - Gas  pressure 

- Defined  as  Apo/  M 

S - Entropy  of  gas 

t - Time  coordinate 

X - Gas  temperature 

a - Gas  velocity  or  projectile  velocity 

U - Internal  energy  of  gas 

V - Specific  volume  of  gas 

X - Position  coordinate 

.i/a 

OC  - Defined  as  q>  p 

P - Entropy  constant  in  semi-empirical  entropic  equation 

- Mean  value  of  specific  heat  in  Noble-Abel  equation 

- Defined  as  I ^EFStcTrvt  (oit  p.fc) 

X - Defined  as  <»yp(p-£) 

p - Density  of  gas 

O’  - Riemann  function,  defined  on  page  4 

(|)  - Entropy  constant  in  semi-empirical  entropic  equation 

The  subscript  "oM  refers  to  initial  conditions  of  the  gae, 
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Fig.  1.  Constants  for  nitrogen  of  the  semi-empirical 
equation  [q)  in  m/sec  (atm)1/*,  / in  Amagat 

p dimensionless]. 


